Higher-order-in-spin interaction Hamiltonians for binary black holes from source 
terms of Kerr geometry in approximate ADM coordinates 
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The Kerr metric outside tlie ergosphere is transformed into Arnowitt-Deser-Misner coordinates 
up to the orders 1/r'* and a' , respectively in radial coordinate r and reduced angular momentum 
variable a, starting from the Kerr solution in quasi-isotropic as well as harmonic coordinates. The 
distributional source terms for the approximate solution are calculated. To leading order in linear 
momenta, higher-order-in-spin interaction Hamiltonians for black-hole binaries are derived. 

PACS numbers: 04.25.-g, 04.25. Nx 



I. INTRODUCTION 



oo 

o 

O 
(N 

>->■ 

. The motion of binary compact objects with proper rotations is an important topic in general relativity. Very likely, 
the first detection of gravitational waves will originate from those sources. Recent achievements of the numerical 
, relativity community in the simulation of merging spinning binary black holes demand deeper analytic work in the 
• treatment of the motion of binary black holes to compare with and hopefully to support the obtained results, P, 

In the analytic treatment of the motion of compact binaries with point-like components the Arnowitt-Deser-Misner 
(ADM) canonical formalism has proved very efficient, 0, 0]; also see Q. It is thus highly desirable to generalize 
that formalism to spinning compact binaries. The aim of this paper is to transfer information from spinning black 
holes at rest, i.e. from Kerr black holes Q, onto the dynamics of binary black holes with spin. For reaching this 
j^jT) goal, at first the Kerr metric is approximately transformed into coordinates which are crucial for the ADM canonical 
formalism. Then the source terms connected with this form of the metric are computed. Finally, new terms for the 
binary interaction Hamiltonian are constructed. 

The paper heavily relies on the diploma thesis by one of the authors, [1]. In the paper, the theory of generalized 
functions is extensively used, @. Throughout the paper, the speed of light c and the Newtonian gravitational constant 
G are put equal to 1. Greek indices will run through 0,1,2,3; Latin indices through 1,2,3. The signature of the metric 
in is +2. 
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^ ■ II. (3+l)-SPLITTING AND ADM FORMALISM 

i> : 

' In this section the ADM formalism is presented as much as it is needed for the purpose of our paper, [loj . 
^ , The (3 -I- l)-decomposition [Tl| of a spacetime manifold gives rise to the following line element: 

^ . ds^ = {N,N' - N^)dt'^ + 2N,dx' dt + -f.jdx' dx^ = -N^dt^ + J^j{dx' + N'dt) (dx^ + N^dt) , (1) 

where N is the lapse function, {Ni = %jN^) the shift vector, and 7.^ the 3-metric of the x'^ = t = const, spacelike 
hypersurfaces. The following relations between these quantities and the metric tensor hold, 



(3) 



with 7'^ being the inverse metric of 7^ , 7*'^7fcj — 5; 
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The metric tensor ([2]) admits a gauge transformation via a diffeomorphism / resulting in a transformation formula 

9U^)--=afi^))^^^ (4) 
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{/"} being the old coordinates depending on the new ones {x^^}- 

The projections of the energy-momentum tensor T^,y related with the (3 + l)-slicing of spacetime read, see, e.g., 

E = nf^n^T^^ = N'^T^^ , (5) 

^ ^l-tn'^T^, = - 1 (To. - T,,N^) = NT^ , (6) 

"^ij = — Tij , (7) 

where = (— A^, 0,0, 0) is the unit vector orthogonal to the t = const, slices and _Lf — (5f denotes the projection 
tensor onto these slices. The source terms in the ADM formalism read [Toj . 

E* = ^E = NV^T^" , (8) 

5"^ V75'' = 7^^V^T0, (9) 

5"^" EE ^N^S'^ = i^^y^V^^Tfei (10) 
and appear in the field equations as follows, 

1 / 1 



^^^TTij — 2""^^ V7-R(3) — IGttE* (energy constraint equation) , (11) 

— TT^'^.j, = SttS'** (momentum constraint equations) , (12) 

dtjij = —'2Nj~^^'^{nij — -TTjij) + Ni-j + Nj;i (evolution equations) , (13) 

dtn^^ ^-N^ ((3)i?,^. _ iy^- i?(3)) + i7V7-'/'7'^ (7r™"7r„,„ - ^tt') - 2N^~'/^i7r"^^7Tj - Itttt^^ ) ^^^^ 
+ v^(iV''^' - 7''^iV'™.„) + (tt*^'^");™ - iV^Tr™^' - N'.^n"'' + 16nS*'' (evolution equations) , 

16^^j^™, EE 7:^^ = ^ (ro, - 7p,r^,7-) 7'f 7^"« (15) 

is the canonical momentum density of the gravitational field, tt = 7y 7r*-' and ";" denotes the 3-dimensional covariant 
derivative. 

The ADM Hamiltonian results from 

-ffADM = 77- / d^Si(7y,j (16) 
lOTT 

after insertion of a solution of the constraint equations subjected to appropriate coordinate conditions. 

A. ADMTT coordinates 

To get a unique solution of the constraint and evolution equations the coordinates have to be fixed by a gauging 
process. The ADM transverse traceless (ADMTT) gauge refers to generalized isotropic coordinates defined by the 
conditions [l3| 
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7y = ( 1 + g'^ ) hi/ , = (17) 
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r TT 



with hfj^ being transversal and tracefree {H, 

be achieved by apphcation of the TT-projection operator Sf^'^'' defined by 



J — , /ifj =0). The TT-part of a (symmetric) tensor of rank two can 



cTTkl 



-{5m- A-^dkdi){5,,- A~^d^dj) . 



(18) 



B. The conformal flatness condition 

Under the conformal flatness condition (CFC) for the 3-metric 7^ , 

7,7^ = ^%-, ^ = l + ^c^, (19) 
which on account of tt*' = results in tt = (maximal slicing), the equations for N, and N'^ read, see, e.g. [T2j . 

A^p = -2n( ij-'E* + ^-^^^i^] , (20) 



A(iVV') =2TrN 



tp-\E* + 2S*) + 



167r 



16n 



(21) 



AN' ^WttNiP-'^S" + 2k'^dj (^-^^ - ^d'dkN'' , (22) 



= 7^ ( S.A^^' + d,N' - -S^.dkN'^ ] , (23) 



where S* = ■^77^5*-' and where Kij is defined by 

±( 

2N V"""' ' 3 
(notice here: KijK^^ — TTy Tr*-'). 

C. CFC+ approximation 

In regard to our Kerr metric calculations, the CFC metric will be generalized to the leading order deviation (LO) 
term from isotropy, 

^CFC+ = + . (24) 

In the ADMTT gauge hfj^ turns out to be the leading order part of hfj^. This hfj^ is a solution of the tensor 
Poisson equation 



Ah^^5™FM, (25) 



where the source F^i is given by 



Fki ~ W.kU. I ~l6Tr^^ (26) 
E* 

with U being the "Newtonian" potential which is a solution of the Poisson equation 

AU^-inE*. (27) 
Considering the CFC equations above for N, tp, and N' only the equation for the lapse function has to be modified: 

A(iVV) = [AiN^)]^rT^, - S^^d^'hJ^^Uu. (28) 
However, this term will turn out be to neglectible for the Kerr metric in our approximation. 
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III. APPROXIMATE ADM COORDINATES FOR KERR METRIC 
A. Transformation from quasi-isotropic to ADMTT coordinates 

We start with the line element in quasi-isotropic coordinates [l^, 



with 



N't' 





(r2 + a2)2 - Aa2 sin^ B ' 



(r2 +a2)2 - Aa2 sin^ 6* ' 

p2/3 ((^2 ^ ^2)2 _ ^^2 gij^2 ^^1/3 

(r2+a2)2- Aa2 sin^ 6* 



and 



= + cos^ ( 



m + a 



r 1 + 



2f 



A = - 2mr + , 
m — a" 



1 



2f 



with (r, 0, 0) denoting the usual Boyer-Lindquist coordinates. Notice that in the limit a 
Schwarzschild metric in spatial isotropic coordinates 



As' 





7n \ 2 




m \ 










2f / 







1^2 , ^2jfl2 ,^2gjj^2 



(dP + f2 jg,2 



(29) 

(30) 
(31) 

(32) 
(33) 

(34) 
(35) 

we get back the 
(36) 



We write the 3-dimensional part of the line element (j29p . = jijdx'^dx^ , in cartesian coordinates via the trans- 
formation 



(f, 9, (p) = \J + y"^ + z2, arccos— , arctan — 
V r X 



^sphere 



which results in 



a;2 + j/2g 2/i ^ ^ 



a;2 + y2 



-dx' 



^2g-2M 



2)2 -t- ?/2 



^dy^ + dz' 



J.2 _|_ y2 



fl-e-2M) dxdy 



(37) 



(38) 



and make a series expansion of (|38p in powers of f = y^E2+2/2 + z2 with the replacement of f by r up to the orders 
1 /r'' and a? . The reason for stopping at these orders is the following: The expansion in ^ corresponds to an expansion 
in powers of ^ with ^ and is clearly nothing else than a post-Newtonian expansion. Notice that it is also possible 
to count in terms of m ^ ^ and a ~ ^ instead of 1/r using m = {M is the mass in Sl-units) and a = 
where the well-known spin S — Mca has been introduced which for rapidly rotating black holes is of the order 1/c. 
So an expansion up to the 4th order in 1/r will suffice to get some of the desired next-to-leading-order (self-) spin 
interaction Hamiltonians for a binary black hole system, 



dl' ^ 



z z z c\ 'A z o222 

m a a x 2a my 2a m y 



2r2 2r3 



dx' 



m \ 4 a'^ 
2^) ^2^ 

2r2 



m^a^ a^t/^ 2a mx 



2a'^m^x^ 



2r3 



8r4 

2 



ma 



2r3 



8r4 



(39) 



dz' 



xya^ 2a^mxy 2a^m'^xy^ . 
2 I '— 1 I dx dy. 
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We will now put this line element into an invariant form with respect to the direction of the spin. For that we 
define 

„ \ / ,x / N a a x ra n , , ^ , ^ 

I , a = (a') = (a,) = - , z ^ — = — — , n = (n') = (n,) (40) 

and use the following relation 

x^dy^ + y^dx^ — 2xy dx dy — [e^ ■ [x x da;))^ 

= -^eijkUpq a" a'' dx'^dx'^ 

= ^ ( aV^ dx^ + 2{a-x){a-dx){x-dx) - {a-xfdx^ '•^"^^ 
a"' \ 

- r'^{a-dxf ~ a^{x-dxy 
with 

,,^,ei,, = 3! <5|;^JJ = + 5\5l5l + 5;5^/f - S^SlS^ - <5;<5^<5f - 5,*5^<5,^ , (42) 

where Cijk is the completely antisymmetric Levi-Civita tensor. Then we find 

1.=!^ +1^+1^+1^. (43) 



where 



lij 




(2) _ 
lij 


1 


(3) 


U 


(4) _ 
lij 





+ (44) 

,^2(a.n)«,.,)^ (45) 

2m?a?nin.i 2m'^aiai 4m2(a-n)a(in,) , , 

+ 4 — ^ + — . ^ ■ (47) 



r3 

7 ^2„2 



8 ^4 ^4 

The index (s) indicates the exact static limit. The transformation of this metric to ADMTT coordinates up to 
the orders l/r* and will be achieved by means of the transformation formula (jASP (the index qiso stands for 
quasi-isotropic) 



7^"^'' = ir + e.3+iT^\^+ ^ti (^s) 



with vector 



^t, 1 a^Ti'' 1 (a-n)a'^ 7 m^a^n^ 7 rn^ia-nYn^ , , 

f — L_\ L 1 !^ L {M)\ 

^ 4 r 2 r 16 r3 ^ 4 r3 ' ^ ' 



This leads to 



itf" - ((l + S)' + 4^ - 5=^) + '■7 + O (<,^ 1 ) (50) 



with 



rprp 7 m?a? m?{a-n)'^ rn?'o?ninj 7v?{a-n)'^ninj 7 iv?aiaj 

^ ~2~^ ^ — 7^ — + ' — 21 '^2~~7 — ■ ^ ' 

The property of leading order hj^^ of being quadratic in is nicely consistent with . In terms of a post-Newtonian 
expansion, we would call hj,^ being of 4th post-Newtonian order, 1/c^. 
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Now the calculations of the lapse and shift functions and the field momentum density are going to be performed. The 
part of the line clement which contains the lapse and shift functions reads (|29p . 



r t 

= {N,N' - N^) dt^ + 2N,dx'dt . 
Using 



(52) 



f2 sin 9 a af"^ sin 9 

the shift vector in quasi-isotropic coordinates is obtained in the form 

, 2 p-2p 

a 



-9 ^ 2 /, ~ • ^ = ^. 2 ^ ^iska'x^dx'' , (53) 



2m 2m 1 (3m^ + ma^) 2m{a-nf- 



f3 2 f4 + -4 j " 



(54) 



together with the lapse function 

m 1 m^ 1 (m^ + ma^) m(a-n)'^ 1 m^ 9 m^a^ 4m^(a-n)'^ 



(55) 



1 



r 2 f2 4 f3 r3 8 f 4 4 



The ADM form of the shift hmction is obtained with the help of the transformation with the vector taken from 
(SHI) 

9^''=9r+9re',.+9ti^'- (56) 

The result reads, again substituting r for f, 

ADM _ .MDM _ f 2m 2m2 3 m3 5m{a-nf ~ma^ \ ■ 



/ 2m 6m^ 21 m^ ^mla-nY - ma?\ , 
^ADM = T + ^ - ^-r + ^ -A ^''^a,nk, (58) 



where use has been made of 

r 2 r 

Next we determine the lapse function in ADM coordinates by means of the very easy formula (|A3p : 



yiDM=(^- — + l4)s^'+0[i,] . (59) 



ADM qiso . qiso ^fe /'^;n^ 

9oo =9oo +9oo,k^ (60) 

with 

qiso 2m 2m^ 3 m'^ Ima? 2m{a-n)^ rn^ Qm?{a-n)^ — rn^a? 

9qo =^'^ + — ^+2^3" +2^3 r3 ~^ ' ^ ' 

resulting in 

ArADM 1 m^ 1 m'^ 1 m** 1 (3m(a-n)'^ — ma^) 1 (5m^a^ — 9m^(a-n)^) 

= 1-7+272"" 4~^8~^2 +2 7^ • (^ ^ 
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Hereof we calculate 



and 



Finally, we calculate tt*-' via 



— ^ 1 2m 5 TO^ / 1 



1 2m 2vn? 3 m^ m^ 3m(a-n)^ — ma^ 2m?a? 

- 1 - — + ^ 4 + — ^ 



with 



6m 3™^ 35m(a-n)^ — 5ma^ 



^3 2 r'^ 



10m{a-n) 



<^ki(pnq)a n g ^ki(paq)a n 



giving us the result 



6m 12m? 15m"^ 35m(a-n)^ — hraa? 



+ 



lOm(a-n) 



(63) 



(64) 



(65) 



(66) 



(67) 



We easily see that the gauge condition on tt'* is already satisfied in this order so there is no need to force to be 
unequal to zero. 

B. Transformation from harmonic to ADMTT coordinates 

First we need the Kerr-metric in harmonic coordinates at least up to the desired order. To achieve this goal let us 
start with Boyer-Lindquist coordinates 



ds^ 



4 [dt - a sin^ ed<ff + ^ \{r^ + a^)d4> - a dil ' + ^dr^ + d^^ 



with 



A = — 2mr + and = r'^ + a? cos^ I 

Rewritten, the line element reads 



, a? SIT? 6b " ^ , ■) 2 sin^ Ob i ^ ,9 9n\ , , 

' —^{/^a-a[rl+a'))d<PB'itB 



A 



^ dr| + p'del + [(r| + a^) - Aa^ sir? 9b] d^l 



(68) 



(69) 



with {tB, tb, Ob, (('b) labeling the Boyer-Lindquist coordinates . 

The harmonic coordinates condition □•xj^ = D^D^^x'^ = (label H for harmonic) gives rise to the equations 

tH=tB, (70) 

and 



XH+iyn = (rs — m + ia)e*'^ sin ( 



In 



rB - r+ 
rB - r_ 



(be aware ip 7^ (fin) 



(71) 
(72) 
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where 

r± = mzL \J w? — a?- (73) 

r\ = XH+yh + z% = {tb - mf + sin^ Ob with (74) 

^'h = r% cosjj = {tb - my cos^ 6b ■ (75) 

These equations are solved perturbatively so that we get the Boyer-Lindquist coordinates depending on the harmonic 
ones, 



a^sin^^H , a^{3 + 5cos{2eH))sm'eH , „ / g 1 \ 



o^smeHCoaOH , 3a'^sin6'i/ cos6'ff cos(26'/f) , ^ / e 1 \ 

^^=^^ — + 8^1 ^ r'^J ' ^ ^ 

= </.^. + + O f a^ -L) . (78) 



2,r]j V 

The line element in harmonic coordinates then reads 

ds^ = goo dt^ + 2go^dt d(f> + dll + dZ^2 , from here on: (r^ ,6HAH)={r,6,4>) (79) 
with the expressions 

. 2m ma^ - 4m^ + Sa^m cos(26') 2mV - 2m^ + 4m2a2 cos(26') 
gQQ = -l + — — + , (80) 



r r 



2am sin^ 9 2am^ sir? 9 am sir? 9 (3o^ - 4to^ + 5a^ cos(26i)) am^ sin^ 6 (3a^ - 2m^ + 5a^ cos(2^)) 
904> = + ;:2 + 2^ ^ ' (81) 

,,9 2m 2m^ 2m^ 2m'^\ , , , so . x9 . o ^ , ,9 2am^sin^0 2m\ , 

= (1 + V + ~ + ~ + — j dr2 + (r + m)W + (r + m)2sm2^d.^2 (^1 + — j drd<^, (82) 

,,9 fa?m,^ .0^ a^m^ „ , 9 a^?™ /l + 3cos(2^) mcos(20)\ ,^9 
d^a^ = ^(3sin^ ^ - 2) - ^(1 + 3cos(20)) dr^ - — + d^^ 

f83^ 

a^?^ . 2. „ to(cos(26») - 3) \ , ,9 2rnVsin(26i) , 

+ -— sin^ ^ -(1 + 3cos(2^)) + — ^ ^-^ ddf ^ , ^ ' drdO . 

2r \ r J 2^ 

dl^ and dZ^2 arc the pure spatial parts of the line element, linear and quadratic in a respectively. Now the harmonic 
3-metric "fij = 7!^"'' + ^ related respectively to dl^ and dl'^2 will be calculated. 
The transformation formula of the metric tensor gives for 7^"^ the result 

(a) 2m m^\ ^ (m^ 2m^ 2m^\ (2m^ 4m3\ u , 

rij =[^ + — + ^)Sii+\^-^ + -^ + ^j n^nj - + j e«(,n,)a'=n' , (84) 

and for ^ 

(a^) _ / '^^^ ~ "^"^ (a-n)^ (a-n)^ - 2m^a^ \ 4m^ (a-n) a(,nj) Sm'^UiUj 

Zm^a^niUj 3to^ (a-n)^ njnj 
From the g^^ part we get the shift vector in the harmonic gauge 



(85) 



,,1, I 2m 2m 2m'^ 5m (a-n) — ma \ ^ u 

N^ =1-^ + — - — + -^-^4 I ^m^^'^' ■ (86) 
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The lapse function in harmonic coordinates turns out to be 



The 3- metric will be brought into ADM- form by means of the transformation formula (jASj): 

7l}°" = 1^J + 7y-, k e + 7fe. e', , + 7fe, + Iki e, (88) 

with vector 

^1 1 , Im? J,,,- 1 ^ > 1 r77?(a-n)a' 3 , , /nr>^ 



The transformation of g^^ is provided by the formula (jA3P : 

SoT^'^sSj+So^.f (90) 

with 

TT 2m 2m? 2m^ 2m^ ma? — Zmla-nY d>m?(a-n? — 2m'^a? , , 

= -1 + — - — + — - — + T ■ (91) 



This results in 



AriM 2m 2m? "im? ma? — "imia-n? '^,m?(a-ri? — 2m^c? , , 

ffo^DM ^ _ 1 ^ _ _ ^ }_ _ ^ ) ^ (92) 



,ADM 



The transformation of is provided by the formula (jA4l 



ADM „H , „H « , „H a (a'W 



which results in 

ADM f 2m 2to^ 3 m,^ 5m(a-n)^ - ma^ \ 
dor"" = [-^ + — - 2— + ^4 j (94) 

giving the same tt*-' as above. 

IV. SOURCE TERMS FOR KERR METRIC IN ADMTT COORDINATES 

The source of the Hamiltonian constraint for a single Kerr black hole presented perturbatively in ADMTT coordi- 
nates is given by 



-A(j) = 



1 , 1 ,2 1 
1 - -6-\ 6^ ( 

8^ 64^ 512 



16TTmS{x) + ^-i + {a-dfl6TTmS{x) + {n'^f + O (95) 



with (a- 9)^ :— a'a^didj — (a-d)^ — ^a?A being tracefree to get rid of the delta- function which appears by calculating 
didj^ and end up with the quadrupole moment of a Kerr black hole (with the specific factor 1/2), which describes 

its deformation. This source is sufficient to reproduce the isotropic part of 7,^ . Notice that (ti'*"')^ is given by the 
Bowen-York expression [l^ and that (j)^5(x) with n G N \ {0} is regularized to zero for a single black hole but will 
contribute for a binary system. From Eq. (95) wc readily read off, 



E* =m[l- -{a-df ]5(x). (96) 



The source for momentum constraint can be calculated to give 

<^ = -/'''akdi {l - l&T:m5{x) + {-]^ + 0,,7r'^ + O (^1) (97) 
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which can be checked by using the formula 



with 



(98) 



(99) 



From the Eq. (97) we find the source expression 



5* 



(100) 



The source for K^'^ in the leading order comes by the Ricci equation Rij = which reads, to leading order, 



Ml. 



TT 



with 



-A"^ IGnm ^1 - ^{a-df ] 6{x) 
4m 2ma^ — 6m{a-n)^ 



TT 



(101) 



(102) 



and 



as well as 



N 



4^ 32^ 256^ 2048^ Wj 



(103) 



(104) 



Notice that our hj, is also the result of the equation 



^^TT _ gTTki f 7(3m^afca; - m'^a'^6ki) \ _ ^TTfe/ 21m^OfcO; 



(105) 



i.e. akUi/r^ is an irreducible part of a metric of the form (ICip to be TT-projected. This result implies that 
the regularized source term in Eq. (26), reading E*{S*k/E*){S*i/E*), vanishes, i.e. [{Sl/E*){SI /E*)]reg = 
[{S'^/E*)]reg[{Si /E*)]reg = 0. Technically one may argue that in the regularization procedure one is allowed to 
apply the "tweedling of products" by Infeld and Plebanski, fl9l|, and to put to zero each term seperately. Such a 
regularization property is particularly valid in dimensional regularization, On the other side, an expansion like 
Eq. (26) would never be expected to represent highly nonlinear terms correctly. So the vanishing of the expression in 
question is quite satisfying. The outcome of our considerations is the vanishing of S*^^ in the given approximation. 



(106) 



A. Leading-order sources for the constraints for boosted Kerr 



To obtain leadingorder source terms with non- vanishing linear momentum we substitute in the constraint equations 
for point particles :5| pS{x) by {p— ^ma x V)(5(a;) according to the a remark in [20] (in the following equation point- 
particle expressions are denoted by the upper index (PP)): 
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-A</. = -Acl^l'Z^ + l^'^P^o^ ) ^'•(^ ^ IGTrSix), (107) 

TT*-^. = -87rp*(5(a;) + 47rTO (a x V)' (5(a;) (108) 

[p = (Pi) a-nd (a x V)* = e^^^ajdk = (a x V)j]. Note that the term in Eq. (107) containing the factor 5/16 is already 
needed at leading order because the integral for the Hamiltonian of the term with factor 1/2 vanishes; however, the 
latter term contributes by iteration via linear momentum independent expressions. 

For a boosted Kerr metric in ADMTT coordinates we can then state the following source expression, to linear order 
in the momentum p and to 6th order in 1/r of the expansion of the unboosted Kerr metric, see Eqs. (95) and (97), 



-A0 = 



1 1 , 1 „ 



1 5 

1 (/) \ p-{ax V) 167r5(a; 



16TTm6(x) + (-1 + ^(j)] (a-dfl6nm6(x) + (n'^f 

V 2 16 / , , 

^ ^ (109) 



6V 512' 
5 

2 ■ 16^^ 

TT^ =ie^"afe9, (^1 - ^(^') 16nmS{x) + (-^ + c^.jT^'' ~ 8np'S{x). (110) 

Here we had to take into account that the substitution in the momentum constraint leads to no new information 
regarding the spin because the expansion of the Kerr metric already gives us the correct linear-in-spin term which 
has no direct connection with the linear momentum. This can be regarded as consistency check for our specific 
substitution. 

B. Source terms in the CFC reduced Kerr metric 

It is straightforward to calculate the source terms in the Eqs. (20)- (22) from our known form of the Kerr metric. 
The outcome reads 

E* ^mfl-^ic^Adix), (111) 



and 



5*=0, (112) 



S*' = -'^e'^'^akdi (^1 - ^(^') Six) . (113) 
All source terms are in full agreement with our previous findings. 

V. SPIN INTERACTION HAMILTONIANS 

The Hamiltonian and momentum constraint equations for point masses are given in the form, see, e.g. Q 



= ^67^Y,{Y'Pa^PaJ+mlf'Sa, (114) 

a 

3 = ^^Y.^"'PaA- (115) 
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The vectors Xa = (x^) (E R denote the position of the a-th point mass, also we define ra '■= x — Xa, ra := \fa\, 
fia '■= fa/ra, and for a ^ b, Tab '■= Xa — Xb- The hnear momentum vector is denoted by = (pai) and the short-cut 
Sa has the meaning S{x — Xa)- The ADM Hamiltonian takes the form 



H 



''^ai Hal 'Si ' 'I j-ji 



•J-O) Pa' ^ij ' T^TT 



(116) 



A. Linear in G spin-spin interaction Hamiltonians 

Our source for the Hamiltonian constraint is generahzed to two black holes by superposition 

- A(/.i = 16^ - i(arar)2^ mi<5i + IGtt - ^{a^Y^ m-2.h ■ 
From this expression the Hamiltonian results in the form 



x = m\-\- mi 



being the rest-mass energy. 

For next order Hamiltonian we must calculate 0i 



(117) 



(118) 



.. = 4(l-l(aT^)2)^+4(l-i(aT^)2)^ 



(119) 



The next order source takes then the form 



-A(/.ii, 



=167r 
+ IGtt 



after Hadamard regularization which disposes of the self- terms like and ^didj6i. 
Integration then gives 



(120) 



Hi, 



1- -miTO2(a2-a2)2 h -^mimiiai-dxY — 

r-12 2 ri2 2 ri2 

-mim2(ar9i)2(a2-92)2 — 

4 ri2 



= -f? jv + Ho 



He 



(121) 



with 



Hn 

Hsl + Hsi 
Hslsl 



77117712 

r-12 

1 mi77l2 

2^ 



[ai-mif -I- 3(a2-ni2)^ - a? - 02) > 



(^(ai-a2) -a^a2j+——^ — \a-^{a2-ni2) + 02 (ai-ni2) j 



3 mir772 

15mim2 (ai-Tii2) (a2-ni2) / 7 
H s I ai-a2 - - (ai-ni2) (a2-ni2) 



r-5 



(122) 
(123) 
(124) 
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Now wc come to the contribution of (tt*^ )^ leading to some new Hamiltonians. The momentum of the binary system 
in the leading order is given by linear superposition 

(7r^^-)2 = (7ri^)2 + (4^)2+27rr7r^^ (125) 

with 

T,^ = -mie'P'diaipdn (^1 - ^{^^y^ 1 - niie^P' dlaipdn (^1 - ^{a^)^^ , (126) 
ni^=nl\1^2). (127) 
The interaction terms will only appear in the object 

27rl^'4^ = 4mim2e'P^aipdudij (^1 - ^{0^^^ ^^e'"' a2qd2sd2j (l - ^{c^y^ ^ 
+ Amim2e'P'aipdudij (^1 - l{a^y^ ^e^'^aa^^Sa^. (^1 - ^{a^y^ ^ 
leading to a Hamiltonian 



= -imim2e'P'aipi9ij5y ^1 - ^{ai-di^^ ("*'' a2qd2sd2j - ^(a2-92)2^ ri2 
- \mxm2e'^^axpdndxj {l - ^{c^^f^ e^"' a2gd2sd2i (l - ^{c^^^ ri2 

= -lm,m2e'P^aipdud,j (^1 - ^{^^^^ e'l' a2gd2sd2j (^1 - ^{a^r^ n2 



(notice du = -821)- 
We then get 



' 12 ^12 



r-12 



_ 35 (ai-ni2)(a2-ni2) 
2 r?2 
= i^5.s?(l--2). 

We sum up expressions 



(128) 



(129) 



HS1S2 = ~rnim2 {ai-a2{di-d2f - {ai-d2){a2-di){di-d2)) ru 

aia2 (ai-ni2)(a2-ni2) ^^^^^ 

— mim2 ^ 6171-11712 3 , 



-f^SiS? = T^'rnim2 {ai-a2){di-d2f{a2-d2Y - {ai-d2){a2-d]){di-d2){a2-d2Y 

= -J^miin2{ai-d2){a2-di){di-d2){a2-d2)'^ri2 

( 3(01-02)0^ ^ 15 (a2-ni2j-(ai-a2) ^ 15 a^(ai ■ni2 j(a2 ■rii2) /-,o1^ 

= mim2 -5 + - -5 + - -5 (i^i) 

\ ^ '12 ^ '12 ^ '12 



14 



1 mi mo 

Z /in 



1^3 (ai-ni2)^ + 3 (a2-ni2)^ - - a2 + 2ai-a2 - 6 (ai-ni2) (a2-ni2)^ , (132) 



which is the same expression as in e.g., [20|, |2]|; also see |22l |23||. 



B. Linear in G and p spin-orbit interaction Hamiltonians 

Taking the partly boosted sources for Hamiltonian and momentum constraint into consideration we will be able 
to compute the spin-orbit Hamiltonian with some higher spin corrections. First we look at the binary source of the 
Hamiltonian constraint 



(133) 



with 



-A<^|^^ i67r ( ( 1 - 1(0^)2') mi + ipi-(ai x d{)\ 5i + 16n Ul - ^{c^)A + x 82) ) 82 



(134) 
(135) 



Integration over — A^j^'' will not contribute to spin-orbit-interaction, but — A0p^ will. So we now have to solve the 
momentum constraint. 



(136) 



with 



47re''='aifc9u 1 - -{ai-diY mi5i - ^T:p\5i = 7r^^,(2 ^ 1) 



(137) 



which gives 



-Pik (Sijdik — - 2Sikdij — ~ 2Sjkdii — + ^dudijdikri] = 712'' (2 ^ 1) 



(138) 



resulting in an interaction contribution (indicated by ~): 



27r;^7r^^ ~ 2 



le'P'^ d{aipdii (l ~ i(ai-5i)2 ) —p2k ( -'26ikd2j— - 2Sjkd2i— + \d2id2jd2kr2 



V2 2- 

le'P'dlaipdii ( 1 - ^{ai-di)'^ ) —p2k { -2Sikd2j— - 25jkd2i— + \d2id2jd2kr2 



f2 



ra 2 



m2e*^'a^a2p92/ ( 1 - 7(a2-92)2 ) — pife ( -2(5iA;9ij— - 2Sjkdu— + ^dudijdikri 



Vi 2- 

m2e-'^''92a2p92; ( 1 - ^(02 •52)^ ) ( -25tkdij^ - 2Sjkdu^ + ^dudijOikVi 



r2 



'n 2' 



(139) 



and leading to a spin-orbit Hamiltonian 
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7^ f {Tr'^Yd' 

mn J ^ ' 



mie'P''p2iaipdi- 82811 ( 1 - j{ai-diy 1 ri2 - m2£'''"piia2pdi- 82821 ( 1 - ^(a2-c?2)2 1 ru (140) 



6^ 



6 



^.•(-.3 X P2) ( - ^ + ^^^) + S2-in.. X p,) (-1 + ^ 

'12 '12 '12 / \'l2 '12 '12 



For some further contributions we calculate 



1 , - — m-1 1 , „ s 1 



1 - 7;{a2-d2)^ — + -P2 02 X ^2 - 
2 J r2 2 r2 



(141) 



which leads to the next order source to be integrated 



-A0II = 167r 



167r 



1^(1)^2 (l-i(a2-92)2j -^,^«P2-(a2Xc?2) 



^P2-(a2 X 82) ^ ^-1 + ^(ai-5i)2j mi5i + ^Pi-(ai x ^i) ^ ^-1 + i(a2-92)2j ^2(52 
+ 16^^Pi-(ai xai)5i |-l + i(a2.02)2j ^ + 167r^p2.(a2 x c?2)<52 | -1 + i(ai.9i)2 J ^ , 



giving rise to a Hamiltonian 



(142) 



ip2-(«2 X 82) y-l + i(ai-9i)2j ^ + lpi.(ai X c?i) ^-1 + \{a2-82Yj ^ 
+ \p,-{a, X (^-1 + i(a2-a2)2^ ^ + ^P2-(a2 X 82) (^-l + i(ai-ai)2^ ^ 
^P2-(«2 X 82) i-l + \{a,-8if\ ^ + \p,-{a, X 8,) (-1 + i(a2-92)2 J ^ 



3 m2 

2 mi 
3 mi 
2 m2 



Pi-(5'ixni2) 1 /3a^Pi-(S'i X ni2) 6(a2-ni2)pi-(5i x 02) 15(a2 •7112)^ Pi • (^i x 7212) ' 



12 



12 



' 12 



' 12 



P2-(52Xni2) 1 /3a?P2-(S'2 X ni2) 6(ai -7112) P2 • (S'2 x ai) 15(ai •ni2)^ P2 • (^'2 x ni2) ' 



rt2 



(143) 



Finally the extended spin-orbit Hamiltonian reads, cf., e.g. [211. 122l|. 
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H 



(a) 
SO 



H 



(b) _ 
SO — 



'->! ['''12 ^ P2) I 


2 


rf2^ 


3 TO2 / 3a2 5'i 


■("•12 


XPl) ^ 


4 mi \ 




^ 


+ 5'2-(ni2 X 




an 

'^12 


3 mi /3af52 


■("12 


X P2) ^ 


4 m,2 V 


'12 





5(ai-ni2)^ 



12 



3to2 e ("12 X Pi) 

o -^l 5 

2 mi 



' 12 



6(a2-ni2) 5r(a2 x Pi) 15(a2-ni2)2 S'r(ni2 x p^) 



' 12 



' 12 



5(a2-ni2) 



' 12 



3 mi 
2 m2 



(ni2 X P2) 



(144) 



' 12 



(ai-ni2) S2-(ai x P2) 15(ai •ni2)2 S'2- (ni2 x P2) 



' 12 



' 12 



We state that the total angular momentum J = L + Si + 82, with spin vectors Sa — rnaa (fl — 1,2) and orbital 
angular momentum L = Xa ^ Pai is conserved in time, 



— = {J,F} = 0, (145) 

where H = Hs,p^ + Hs,p, + Hs^p, + Hs^p^ + iJsiSs + Hg2 + iJ^s + Hg3p^ + i/^ipj + i?sf S2P2 + -^siSipi + ^sf si + 
Hg^gz + Hg^gs, by making use of the standard Poisson brackets 

{x^, pbj} = SijSab , {^ai, S'bj} = iijkSakSab , and zero otherwise. (146) 

Also conserved are the absolute values of the spin vectors, i.e. = const. This is consistent with results from 
the literature [l^, [1^, where it is shown that the tidal friction enters on the spin-interaction-separation scaling of 
(l/''i2)^ only. It is also consistent with Ref. '26*1 where tidal deformations for black holes are identified to occur at 
5th post-Newtonian order for the first time, whereas our developments do not exceed the 4th post-Newtonian level. 
Obviously, our gravitating mass parameters nia do not change in relation to the irreducible mass parameters Ma as 
given by ml = Ml + Sl/AMl, see 



VI. CONCLUSIONS 



In this paper several new interaction Hamiltonians for spinning binary black holes have been derived, Hg2 g_^p_^ + 
HsIpi + ^sfp2 + ^SfS2P2^ Hg2g2, H g^g3 -\- Hg^ga. Thcsc Hamiltoiiians generalize the previously known ones, 
Hsipi + Hg^p^ + Hg^p^ -\- Hg^p^, Hg^g^, Hgj + Hg2, through further nonlinear terms in the spins. If the spins 
are counted of the order {1/cf, the new Hamiltonains are of the order 1/c^ higher compared with the previ- 
ously known ones. Other Hamiltonians which are of the same higher order in 1/c, i.e. 1/c^, read Hg-^p.^(^p2_\_Q^ + 

Hs2Pl{p'^+G)l ■H^5i52(p2 + G), Hg2^^p2j^Q) + Hg2(^p2j^Q), H g2 g^p^ + H g2 g^p^ , Hp^g3 + Hp^g3, Hgf + H gi ThcHamil- 

tonian Hg^p^(^p2^Q-^ + Hg^p^(^p2^Q-^ has recently been found by Damour, Jaranowski, and Schafer [2^, the Hamil- 
tonian Hg^g^(j,2^Q^ may possibly result from a recent paper by Porto and Rothstein, [2^. The Hamiltonians 
Hsl(p'+G) + ^s|(p2+G), Hg2g^p^ + Hg2g^p^, Hp^g3 + Hp^gs and Hg4 + H gi are still unknown. In forthcoming 
papers, expressions for all the mentioned Hamiltonians are expected to be presented as well as applications of them 
within the context of gravitational wave astronomy. 
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APPENDIX A: COORDINATE TRANSFORMATION 

We start with an infinitesimal coordinate transformation with translation vectors e.g. [si 
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x''ix"^)=x"'+Cix'^), (Al) 

whereby a;" are the old coordinates and x'" are the new ones the metric tensor is expressed in. The general transfor- 
mation formula together with Taylor expansion leaves us with an expanded transformation formula 

dx" dx^ 

= 9a0{xn{S"^+C^,)iS^.+f..) (A2) 
= {g^x'^ + gap,xix'ne + ■■■) {S'^.S^. + , + -5^^^^, . + , .) 

Under the assumption that the metric tensor is stationary we specialize this formula for the components: 

ffoo = 500 +500,,: f + . . . (A3) 

9m = 9m+9m.]^^ +gaoC^ + ■■■ (A4) 

+9o^e.+90Je,^+gme,^e.+9oke,^e.+9ake,^e. + ■• • (as) 



APPENDIX B: USEFUL FORMULAS 

From the Ref. [3l| we take the formula 

A(r^nL) (A-;)(A + /+l)r^-2nL (VAeC!) (Bl) 

with radial distance r and the short-cut notation fiL nf^'TTThi^, whereby — x'^/r is the normal unit vector and 
the hat indicates tracefreeness according to formulas 

"■i = X! ^^^'^ (21 - • • ■ '5i2fc-l»2fc"»2fc + l...*l} ' (B2) 

fe=0 ^ '" 

_ ^ ^ ^fc (2/-4fc-H)!! ^ A - 

'^i ^ 2^ fO/ — 2k + IV' ' ' ' *2fc-l»2fc"-i2fc + l --il} 1 (^3) 

k=0 ^ 

where A^^-^ i^y is the non-weighted sum J2aes ^iai---icTi with S as the smallest set of permutations (1 . . . /) that fully 
symmetrizes Afj J , J -J in Ji . . . ii. 

Inverting (|B1[) leads to the very useful formula 

APPENDIX C: TT-PART OF KERR 3-METRIC IN THE LEADING ORDER 

We make an ansatz for the 3-metric that is to be TT-gauged (6, c,d,e, f,9 being constants) 

= ^^^^0 + c + d + e ^ + 2/ + 9^:^ ■ (CI) 

An infinitesimal coordinate transformation is given by 

lU=1^3+i^,3+^3,^ (C2) 
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with being set as 



f = A- 



B 



^ r hC- ^ ' 



pa — 1 



where A, B, C denote constants. 

The transformed metric then becomes 



jl^^ib + 2A)-—S., + {c + 2B) 



m?{a-ny 



(5y + {d- 2aA) 



(e-2(a + 2)B) 



m'^{a-n)'^ninj 



(2/ + 4B - 2aC) 



m?aiaj 



+ (3 + 2C)— ^ 



The coefHcients will be chosen so that 7- gets TT-gauged. Therefore, we use the condition 



The coefRcients then result in 



A 
B 
C 



For a = A we find: 



•TT 



2d{a - 5)(a - 3)(a + 2) - (a ~ 6)(2e + (2 + a)(2/ + go)) 
4(a-5)(a-3)a(a + 2) ' 
(2 + a)(2/ + ga) + 2e(14 + a{a ~ 8)) 
4(a-5)(a-3)(a + 2) 
2e(6 + (a - 6)a) + (2 + a)(ffa(2a - 9) + 2/(6 + (a - 6)a)) 
2(a-5)(a-3)a(2 + Q;) 



rr?(a-ri)^ m^a^riinj m'^{a-n)'^ninj rn^aiaj 

-^Sij + 2 -g (5y + 2 -5 6 -J ^ -J 



with 



(C3) 



(C4) 



(C5) 



(C6) 
(C7) 
(C8) 



(C9) 



77 = e + 6/ + 12g. 



(CIO) 



APPENDIX D: METHODS OF REGULARIZATIONS 



1. Partie finie regularization according to Hadamard 

Let us consider / being a real function defined in an environment of the point Xq € TR^, except in this point where 
/ is singular. We define a family of complex-valued functions fn as follows: 



U : €3e ^ U{e):= f{xQ+en) e C. 
We expand /„ in a Laurent-series around e = 0: 



(Dl) 



(D2) 



The regularized value of / at Xq is defined as the coefficient at e'^ in the expansion (|D2p mean-valued over all unit 
vectors n, [HIMIll, 



freg{xo) := ^ ^ df^ao(n) 



(D3) 



This formula can be used to calculate integrals with delta-distributions. We define 

J d^X f{x)S{x - Xa) fregiXa) , 

which provides us with a formula for calculating Poisson integrals of the form 

\ a ) \ a ) a a ^ 

Two examples shall clarify this. 

\r I in \r / r 



because 



/=- /„ = /(0 + en) = i ^ (- 

r e \r 



x=0 



(2) A-^ (^^d,d,S{x)^ = . 

To prove this relation we make use of the identity 



-didj6{x) — didj [-5{x) \ ~ S(x)didj- — diS(x)dj- ~ djS(x)di- , 
r V'"/ r r 

d.,S{x)d.j^ = di (^6{x)dj^^ ~ S{x)didj^ , 

^didjS{x) = didj ^i(5(a;)^ + d{x)didj^ - di (^S{x)dj^^ - dj (^d{x)di^^ , 



and conclude 



( ^d,dj6{x) \ = d.djA-^ (~^(^)) +^"^ U(cc)a,9ji j - d,A'^ Ui^)^]^. ) -c>o^'^ ( Kx)di 



47r V / x=o r 47r Vr^ / x=o r 47r V r*^ / 



because of absence of an e^-coefficient. 



2. Riesz's formula 



From the Ref. [35j we adapt the famous formula 



reg 



e.g.. 



3,/mr(^)r(-^ 
r(-f)r(-|)r(^) 

7d3.A" 



Q+/9+3 
12 I 
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